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ABSTRACT 

We present results of the study of persistent high mass X-ray binaries (HMXBs) in the 
Milky Way, obtained from the deep INTEGRAL Galactic plane survey. This survey 
provides us a new insight into the population of high mass X-ray binaries because 
almost half of the whole sample consists of sources discovered with INTEGRAL. It is 
demonstrated for the first time that the majority of persistent HMXBs have super- 
giant companions and their luminosity function steepens somewhere around ~ 2x 10 36 
erg s _1 . We show that the spatial density distribution of HMXBs correlates well with 
the star formation rate distribution in the Galaxy. The vertical distribution of HMXBs 
has a scale-height h ~ 85 pc, that is somewhat larger than the distribution of young 
stars in the Galaxy. We propose a simple toy model, which adequately describes gen- 
eral properties of HMXBs in which neutron stars accrete a matter from the wind of 
the its companion (wind- fed NS-HMXBs population). Using the elaborated model we 
argue that a flaring activity of so-called supergiant fast X-ray transients, the recently 
recognized sub-sample of HMXBs, is likely related with the magnetic arrest of their 
accretion. 
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1 INTRODUCTION 

A formation and evolution of stars and binary systems are 
very slow processes (typical timescales of > 10 6 yr) and can 
not be directly traced. The only way to understand them 
is to study populations of sources. Properties of an ensem- 
ble of sources at different stages of their evolution do have 
characteristics which can be measured and compared with 
predictions of different models. 

X-ray sources in our Gal a xy w ere discovered quite 
long ago l|Giacconi et alj 1 19621 . Il97ll ) and studied inten- 
sively in different ways during ~ 40 years. Advances in 
a sensitivity and angular resolution of current X-ray tele- 
scopes allow now to detect and co unt such objects in a 



varie ty of outer galaxies (see e. g. Trinchieri fc Fabbianol 



19911: iPrimini. Forman. fc Jones! 1 19931 : iGilfanovl |2004 



Kim fc Fabbianoll2004) 



In spite of this observational progress the properties of 
populations of different types of X-ray sources are not fully 
understood yet. The main reason for this is a scarce informa- 
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tion about physical properties of large samples of HMXBs. 
Objects in outer galaxies can be counted and their X-ray 
appearances measured, but physical properties of detected 
sources (first of all, types and ages of their stellar compo- 
nents, orbital periods, etc.), which can help to understand 
their formation and evolution, rarely can be obtained (typ- 
ically from optical or infrared instruments) due to their ex- 
treme faintness. 

On the other hand - physical properties of stellar bi- 
naries in our Galaxy can often be well measured, but the 
uniform highly sensitive surveys of such objects were absent 
until recently. A strong absorption of soft X-rays (< 2 keV) 
in the interstellar medium or near binaries themselves have 
not allowed to obtain an unbiased sample of s ources in the 
Gala ctic plane from the ROSAT all sky surv ey |Voges et al 
1999h . Surveys, performe d with the ASCA (jSugizaki et al 



20011) and XMM-Newton (|Hands et alj|2004 ) observatories, 
are somewhat free from these limitations, but both observa- 
tories covered only small areas of the Galactic plane. 

X-ray instruments, which operate in the harder X- 
ray energy range, where the interstellar photo-absorption 
does not play an important role, and can cover the 
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whole Galaxy, have typically a poor angular resolution 
and thus are affected by a source confusion in the 
Galactic plane (e.g., UHURU.lpbrma.net al.lll97Sl : RXTE, 
iMarkwardt et al.ll200cl ; iRevnivtsev et al.ll2004 l. Samples of 
sources, collected over all history of X-ray astronomy (e.g., 
iLiu. van Paradiis. fc van den Heuvell 120071 1 can not to be 
taken as a representative for statistical and physical studies 
of populations due to their non-uniformity respect to the 
flux, detection criteria, etc. 

A systematic survey o f the Galaxy in 2003- 2011 with the 
INTEGRAL observatory (| Winkler et alj|2003j ) in the hard 
X-ray energy range (> 17 keV) with the moderate angular 
resolution (~ 12') allowed us for the first time to overcome 
all these difficulties and to perform virtually unbiased search 
for X-ray b inaries in the Milky W ay with an unprecedented 
sensitivity l|Krivonos et al.ll2012t ). 

Previous observations of galactic sources with the IN- 
TEGRAL observatory have proved to be fruitful in a va- 
riety of fields: for explanation of the origi n of the Galac- 
tic R idge X-ray emission in hard X-rays jKrivonos et al.l 
2007): for systematic discoveries of strongly photoabsorbed 
high mass X-ray bi naries and study of their distribution in 
the Galaxy (see e.g. Courvoisier et al. 2003; Lutovinov et al 



INTEGRAL/1SGRI 17-60 keV 
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20051 : IBodaehee et all 120071 : iLutovinov et al l 120071 : IChatv 
20081 : IColeiro fc Chatvll201ll ; IBodaehee et~ai1l2012t ) ; for con- 
firmation of the presence of the br eak in the luminosity func - 
tion of low mass X-ray binaries (jRevnivtsev et al . 2008a), 
which is likely related with changes of the evolutionary type 
of st ellar companions in these systems (|Revnivtsev et all 

Eon]). 

In this paper we present the most sensitive unbiased 
flux-limited sample of non-transient high mass X-ray bina- 
ries in our Galaxy. Our main goal is to establish and under- 
stand general properties of this population. This knowledge 
is a very important for different reasons: for understanding 
of properties of compact objects, their formation and evolu- 
tion; for explanation of the observed behaviour of different 
types of HMXBs (in particular, supergiant fast X-ray tran- 
sients) ; for interpretation of observations of variety of galax- 
ies, which are now became possible with new generation of 
X-ray telescopes (in particular, to calculate a contribution 
of HMXBs to observed luminosities), etc. 



2 INTEGRAL SURVEY OF THE GALACTIC 
PLANE. THE SAMPLE OF HMXBS 

In our study we use result s from the INTEGRA L Galactic 
plane survey, presented in iKrivonos et al.1 (|2012D . The sen- 
sitivity of this survey is typically 10~ n erg s _1 cm~ 2 in 
the 17-60 keV energy band, which ensures the detection of 
sources with luminosities > 10 35 erg s _1 within a half of the 
Galaxy ( < 9 kpc from the Sun) and > 5 x 10 35 erg s _1 over 
the whole Galaxy ( < 20 kpc from the Sun) . 

One of advantages of the INTEGRAL Galactic plane 
survey is the multiple coverage of virtually all galactic 
sources. During the period of operation all parts of the 
Galaxy were observed tens of times, that allow us to quan- 
tify a long term behaviour of a majority of sources. This 
is important because the goal of our work is to understand 
properties of the population of high mass X-ray binaries, 
which are typically registered (or can be registered) in ex- 
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Figure 1. Three upper panels. Examples of light curves of 
three sources of different types: persistent (GX 301-2), super- 
giant fast X-ray transient (IGR J17391-3021) and Be-transicnt 
(4U 0115+63). Bottom panel. Normalized distribution of values 
of fluxes, measured for these three sources in single revolutions 
(with the duration of about 3 days) of the INTEGRAL spacecraft 
(histograms). Dashed lines indicate median fluxes registered from 
sources, dashed-dotted lines - mean flux values. Red color corre- 
sponds to the source GX 301-2, black one - to IGR J17391-3021, 
green one - to 4U 0115+63. 



isting and future surveys both our Galaxy and distant galax- 
ies. Therefore in our work we concentrate only on persistent 
sources. Transiently appearing HMXBs typically have small 
duty cycles and thus they should be a minority among the 
population of HMXBs detected in snapshot observations of 
galaxies. 

As a first step, among sources detected on the average 
map of the Galactic plane we exclude all sources which are 
known to be not high mass X-ray binaries. The remaining 
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Figure 2. Sensitivity of the INTEGRAL survey (at the 4.7a level 
in the 17-60 keV energy band) over the whole Galaxy (solid line). 
Coordinates and fluxes of 26 non-identified persistent sources are 
shown by circles. Uncertainties correspond to la. Dashed line 
denotes the survey flux limit (see text for details). 

objects are either HMXBs or unidentified sources. In to- 
tal, 75 HMXBs with a confirmed nature and 33 sources of 
an unknown origin (including HMXB candidates) at lati- 
tudes \b\ < 5° we re detected during th e INTEGRAL Galac- 
tic plane survey (|Krivonos et alj|2012t l. 

Then we exclude transient sources from our sample. 
There is no firmly established concept of a transiency, there- 
fore this is a non-trivial task. To illustrate it, some examples 
of a different long term behaviour of sources are shown in 
Fig[T] where one can see light curves (upper panels) and 
distributions of fluxes (bottom panel) of three sources, rep- 
resenting three big families of high mass X-ray binaries: 
(quasi-)persistent sources, Be-transients and supergiant fast 
X-ray transients (SFXTs). We quantify the "persistency" 
of a source via a ratio of two flux values: the median flux 
value Fmedian and the mean flux value (F). For our main 
source sample we accept only sources which have the ratio 
FmcdianA-F) > 0.5. This resulted in 54 confirmed HMXBs 
and 26 unidentified sources. It is necessary to note, that sev- 
eral SFXTs or candidates, which are named "trantients" , 
still passed the selection criteria. It can be connected with a 
high duty cycle for some sources (like IGR J16479-4514) or 
due to a luminosity dynamic range not particularly high for 
some sources in the INTEGRAL energy band (this could be 
the case of the so called intermediate SFXTs IGRJ17354- 
3255, IGRJ16418-4532, AX J1845. 0-0433). 

In order to further increase the identification complete- 
ness of our sample we raise the flux limit of the survey, as 
virtually all non- identified sources are faint. The sensitivity 
of the survey is shown as a function of the galactic longi- 
tude in Fig[2] It was calculated as the average value of the 
INTEGRAL /IBIS /ISGRI sensitivity over Galactic latitudes 
—5° < b < 5° for each longitude direction. Circles denote 
fluxes of non-identified sources with their uncertainties. Two 
important facts are clearly seen from this figure: 1) the sur- 
vey sensitivity is not uniform over the Galaxy and is much 
better in its inner part (directions with |Z|<100°); 2) all 
non-identified sources are faint. Thus, if we increase the sur- 
vey flux limit up to 0.7 mCrab (10 -11 erg s _1 cm -2 in the 
17 — 60 keV energy band) over the Galactic longitude range 



— 107° < I < 136° and 1.5 mCrab in the remaining part 
of the Galaxy there will be no non-identified sources in our 
sample. At the same time, the list of known HMXBs will 
shorten to 48 sources. 

This sample contains 37 persistent sources with accret- 
ing neutron stars and known distances, five HMXBs with 
unknown distances and six sources of a peculiar nature 
(HMXBs harboring black holes and 7-loud HMXBs). 

The list of selected persistent HMXBs (and non- 
identified) sources with their main parameters is presented 
in Table Q] its contents are described below. 

Column (1) "Name" - source name. 

Columns (2,3) "I, b" - source galactic coordinates, lon- 
gitude and latitude, respectively. 

Column (4) "Luminosity" - time-averaged source lumi- 
nosity in the 17 — 60 keV energy band. For sources with 
unknown distances the flux in mCrabs in this energy band 
is mentioned. 

Column (5) "Distance" - distance to the source in kpc. 

Column (6) "P or b" - orbital period of the system. 

Column (7) "Class" - optical class of the normal com- 
panion in the binary system. In some cases two possible type 
are indicated. 

Column (8) "References, notes" - references and 
alternative source names. References correspond to the 
distance, class and orbital period measurements. 

The spatial distribution of HMXBs from Table [1] are 
shown in Fig[3] 



3 LUMINOSITY FUNCTION AND SPATIAL 
DISTRIBUTION OF WIND-FED HMXBS 

In our Galaxy there are (and it can be seen from the Ta- 
bl^TJ several types of persistent X-ray binaries with massive 
companions. The most numerous population of them is bi- 
naries, in which the neutron star accrete a matter from the 
stellar wind of the massive companion (wind-fed systems). 
Other types of persistent X-ray binaries with massive stars 
include: 

• binaries in which neutron stars accrete matter due to 
the Roche lobe overfl ow of the giant companion star (likely 
Cen X-3, see e.g. lLamers. van den Heuvel. fc Pettersonl 

Glzi); 

• binaries with accreting black holes (CygX-1, likely 
CygX-3); 

• binaries with a super-Eddington regime of accretion, in 
which the central engine is obscured by the accretion disk 
and only jet is visible in X-rays (e.g., SS 433. lFa,brikall2004l ; 

• binaries where an X-ray (and gamma-ray) emission 
originates as a result of a non-thermal emission of particles 
accelerated in colliding winds of components (PSRB1259- 
63, LS5039, LSI +61 303, Eta Carinae, etc.). 

These classes of sources have only a few representatives 
in our Galaxy, which makes a study of their statistics im- 
possible. Therefore we will concentrate in this paper only on 
binary systems with the wind-fed accreting neutron stars. 

A distribution of HMXBs over their luminosities (lumi- 
nosity function, LF) is the simplest, but still an informative 
characteristic, which can be calculated from the sample of 
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1. List of persistent galactic high-mass X-ray binaries detected by the INTEGRAL observatory 



HMXBs with fluxes > 10 11 erg s 1 cm 2 and known distances 



Name 


I, 


b. 


Lx,l7-60keV> Distance, 


Porb 5 


Class 


References 




deg 


deg 


10 35 erg s 1 


kpc 


days 






Vela X-l 


-96.93 


3.93 


0.868 ± U.003 


1.1 


8.96 


B0.5Ib 


1, 2, 3 


3U 1022-55 


-74.64 


1.49 


n si i _i_ n n q q 
0.311 ± 0.U33 


5.0 




BOIII-Ve 


4, 5 


Cen X-3 


-67.90 


0.33 


o A c: i _l_ n c\A 1 
24.01 ± 0.U41 


5.7 


2.09 


06-7II-III 


6, 7, 8 


IGR J11305-6256 


-66.05 


-1.48 


U.299 ± 0.012 


3.0 




BOIIIe 


9, 10 


IGR J11435-6109 


-65.12 


0.68 


Q 1 UK _l_ n 11(17 

3.165 ± 0.09 / 


8.6 


52.46 


B2III or B0V 


11, 12 


A 1145.1-6141 


-64.50 


-0.02 


oa 1 A _l_ n nncr 
20.14 ± 0.U95 


8.5 


14.4 


B2Iae 


13, 14 


X 1145-619 


-64.38 


-0.24 


0.271 ± 0.012 


3.1 


187.5 


BlVne 


15, 5, 16 


1ES 1210-646 


-61.13 


-2.31 


U.1U5 ± U.U11 


2.8 




B2V 


11 


GX 301-2 


-59.90 


-0.03 


61. Do ± U.Ulu 


3.5 


41.5 


Blla+ 


17, 18, 19 


1RXP J130159.6-635806 


-55.91 


-1.12 


U.7do ± 0.U41 


5.5 




09V or B1III 


20, 9 


4U 1416-62 


-46.98 


-1.57 


n aoq _i_ n n a o 
U.4oo ± U.U45 


6.0 


42.12 


BlVe 


21, 5, 22 


4U 1538-522 


-32.58 


2.16 




4.5 


3.73 


B0.2Ia 


23, 24, 25 


IGR J16207-5129 


-27.54 


-1.05 


1 con _i_ n dAi 
l.Doy ± U.U41 


6.1 


9.73 


07.5 


26, 27, 28 














BHa 


26 


IGR J16195-4945 


-26.44 


0.33 


n ac\a _L n noo 
U.4y4 ± U.Uzz 


4.5 




Blsg 


29, 27 


IGR J16318-4848 


-24.38 


-0.44 


U.o9U ± U.UUz 


1.6 




sgB[e] 


27 


IGR J16320-4751 


-23.67 


0.16 


o nncr _i_ n ni q 
z.yyo ± U.Ulo 


3.5 


8.96 


081 


27, 30 


AX J 163904-4642 


-21.99 


0.07 


i c\a a _L n nc: 

/.yb4 ± u.lzo 


10.6 




BIV-V 


31 


IGR J16418-4532 


-20.81 


0.49 


n a r 7'7 _L n 1 m 

y.oYY ± u.iyi 


13.0 


3.75 


08.5(sg?) 


31, 27, 32 


IGR J16465-4507 


-19.94 


0.13 


1.714 ± 0.100 


9.4 


30.32 


B0.5I 


27, 33 














09.51a 


26 


IGR J16479-4514 


-19.84 


-0.12 


0.445 ± 0.009 


2.8 


3.32 


08.51 


26, 27, 34 














09.5Iab 


26 


IGR J16493-4348 


-18.62 


0.57 


6.819 ± 0.260 


15.0 


6.78 


B0.5Ib 


35, 36, 37 


OAO 1657-415 


-15.63 


0.32 


A Q 7H _1_ n HKQ 

48. 1 1 ± 0.05o 


7.1 


10.4 


B0-6sg 


38, 39, 40 


4U 1700-377 


-12.24 


2.17 


12.98 ± 0.004 


2.12 


3.41 


06.5Iaf+ 


41, 42, 43 


EXO 1722-363 


-8.50 


-0.35 


3.996 ± 0.032 


6.1 


9.74 


B0-B1 la 


27, 44, 45 


AA J 1749. 1-2733 


1.58 


0.06 


3.285 + 0.128 


13.5 




Bl-3 


46 


AX J1749.2-2725 


1.70 


0.11 


2.760 + 0.137 


14.0 




B3 


46 


IGR J18027-2016 


9.43 


1.03 


10.24 + 0.140 


12.4 


4.6 


Bib 


47, 48 


IGR J18214-1318 


17.67 


0.48 


1.377 + 0.076 


8.0 




B0V-O9I 


49 


AX J1845.0-0433 


28.14 


-0.66 


0.255 + 0.015 


3.6 




09.51 


50 


XTE J1855-026 


31.07 


-2.09 


14.44 + 0.114 


10.0 


6.07 


BOIaep 


51, 52, 53 


X 1908+075 


41.89 


-0.81 


9.658 + 0.047 


7.0 


4.4 


07.5-9.5sg 


54, 55 


4U 1907+097 


43.74 


0.47 


4.389 + 0.024 


5.0 


8.38 


08-9Ia 


56, 56, 57 


IGR J19140+0951 


44.29 


-0.46 


1.687 + 0.012 


3.6 


13.56 


B1I 


47, 27, 58 














B0. 51 


26 


SWIFT J2000.6+3210 


68.98 


1.13 


2.108 + 0.092 


8.0 




early BV or mid Bill 


59 


4U 2206+543 


100.60 


-1.10 


0.852 + 0.010 


2.6 


9.57 


09.5V 


60, 61 


1A 0114+650 


125.71 


2.55 


6.642 + 0.063 


7.2 


11.6 


BHa 


62, 63 


RX J0146. 9+6121 


129.52 


-0.80 


0.096 ± 0.010 


2.5 




BlVe 


64, 5 


HMXBs and candidates to HMXBs with fluxes > 10 


erg 


— 1 —9 

s cm 


and unknown distances 




IGR J10100-5655 


-77.76 


-0.67 


0.713 ± 0.086 a 






early giant 


70 


AX J1700. 2-4220 


-16.23 


-0.03 


1.124 + 0.070" 




44.03 


Be 


71 


IGR J17200-3116 


-4.99 


3.34 


1.500 ± 0.048 a 






HMXB 


70 


IGR J17354-3255 


-4.54 


-0.26 


0.958 ± 0.045 a 




8.448 


HMXB 


72, 73 


IGR J17586-2129 


8.01 


1.36 


1.502 ± 0.052 a 






HMXB 


74 






HMXBs with fluxes < 10" 


11 erg s' 


_1 cm -2 






IGR J00370+6122 


121.21 


-1.42 


0.084 ± 0.010 


3.0 


15.665 


BN0.5II-IIIb 


68, 69 


IGR J14331-6112 


-45.15 


-0.76 


1.053 + 0.135 


10.0 




Bill or BV 


59 


IGR J16283-4838 


-24.66 


0.08 


3.032 ± 0.339 


17.6 




OBsg 


65 


IGR J17544-2619 


3.24 


-0.32 


0.130 + 0.010 


3.6 


4.93 


091b 


27, 66, 67 


IGR J18410-0535 


26.78 


-0.23 


0.113 + 0.011 


3.2 




BHb 


26 


IGR J22534+6243 


109.87 


2.88 


0.417 ± 0.072 a 






HMXB 


75 
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Table 1. (continue) 



HMXBs with black holes 



Name 


J. 

deg 


b, 
deg 


Lx,17-60keV > 

10 35 erg s _1 


Distance, 
kpc 


days 


Class 


References 


SS 433 
Cyg X-1 
Cyg X-3 


39.69 
71.34 
79.85 


-2.24 
3.07 
0.70 


3.786 ± 0.030 
38.96 ± 0.005 
101.2 ±0.057 


5.5 
1.86 
7.2 


13.1 
5.6 
0.2 


Asg 

09.7Iab 
WR 


76, 77 
78, 79, 80 
81, 82, 83 


7-loud HMXBs 


PSR B1259-63 
LS 5039 
LSI 61 +303 


-55.82 
16.87 
135.67 


-0.99 
-1.29 
1.07 


0.119 ±0.010 
0.075 ± 0.007 
0.143 ±0.013 


2.3 
2.9 
2.0 


1236.7 
3.906 
26.496 


B2e 

06.5Vf 

BOVe 


84, 85, 86 
87, 88, 89 
90, 91, 92 



flux in the 17-60 keV energy band, in mCrabs 



Refere 
C 
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sources. It can be easily done for outer galaxies, but in the 
case of the Milky Way, while calculating the HMXBs LF, 
we should make an appropriate correction for an incomplete 
coverage of the Galaxy at different X-ray luminosities. In 
other words - we should correct for the fact that sources 
with given luminosities are detectable for the INTEGRAL 
survey only within some distance limits. 

The simplest way to do such a cor rectio n was adopted 
by iGrimm. Gilfanov. fc Sunvaevl (|2002p and IVoss fc Aiellol 
(|20ld ). where authors assumed some particular volume den- 
sity distribution of HMXBs over the Galaxy. In the former 
work it was assumed that the HMXBs volume density dis- 
tribution has a disk shape with certain parameters, in the 
latter paper authors suggested that HMXBs are distributed 
similar to the stellar mass in the Galaxy. 

In this work we want to measure the HMXBs density 
distribution over the Galaxy rather then assume it. There- 
fore we have adopted here a different approach. 



3.1 Volume limited samples 

As a first step we divided the whole available luminosity 
range of HMXBs into two intervals - above the luminosity 
2 x 10 35 erg s _1 (taking into account the adopted flux limit 
of our survey 10 -11 erg s _1 cm -2 this sample is complete up 
to ~ 13 kpc from the Sun) and above the luminosity 2 x 10 34 
erg s _1 (complete up to ~ 4.1 kpc from the Sun). For such 
volume-limited samples it is not needed to make a luminosity 
dependent correction and the luminosity distribution can 
be calculated directly. Results of these calculations - two 
portions of the luminosity function of HMXBs - are shown 
in Fig[4]by histograms. For the faint part of the luminosity 
function we present both observed number of sources in the 
volume limited sample (at distances < 4.1 kpc from the 
Sun, dotted line) and the same histogram multiplied by a 



factor of 6.6, which is a re-normalization, calculated from 
the measured density distribution of HMXBs (see below). 

This plot provides us an indication that the luminosity 
function of HMXBs in the whole luminosity range (10 34 — 
10 37 erg s _1 ) is not following a single power law, but rather 
is curved at luminosities around (0.4 — 2) x 10 36 erg s . The 
maximum likelihood approximation of the LF with simple 
power law functions at L > 4 x 10 35 erg s _1 and L > 2 x 
10 34 erg s _1 gives the best fit values of their slopes 7f a i nt = 
1.49 ± 0.21 (34.3 < logL^ < 36.5) and 7 br i g ht = 2.0 ± 0.3 
(35.6 < logL^, < 36.5). From a purely statistical point of 
view this difference in slopes is ~ 2a significant. 

It is important to emphasize that distances to sources 
(Table[T]) are known with a limited accuracy due to different 
reasons. In general, it is reasonable to assume that they have 
an average accuracy not worse than ~20%. It means that 
the luminosity of HMXBs will be uncertain with a factor 
of ~40% for the known average flux values. To estimate 
a possible influence of these uncertainties on our results we 
have performed simple simulations. We have varied distances 
around values, presented in Table [T] assuming a Gaussian 
distribution with cr=20% of the source distance. Systematic 
uncertainties for the LF slopes calculated by this way are: 
A7f a i n t ~ 0.06, A7bri g ht ~ 0.1, respectively, i.e. well within 
the statistical uncertainties. 



3.2 Luminosity function of the whole sample 

In order to calculate the luminosity function of the whole 
sample of HMXBs we should take into account the fact that 
sources with given luminosities are detectable for the INTE- 
GRAL survey only within some distance limits. In order to 
do this we need to determine their surface density distribu- 
tion. 

Given the limited accuracy of distances to sources in 
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Figure 3. Maps of the inner part of the Galactic plane, obtained with INTEGRAL/IBIS in the 17-60 keV energy band. All persistent 
sources from Table [l] arc marked with circles and names. 
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Table 2. Best fit parameters of the luminosity function of 
HMXBs and their spatial density distribution 
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Figure 4. Luminosity functions of the wind fed accreting HMXBs 
in the Galaxy. Two black solid histograms represent luminosity 
functions within volume limited samples (see text for details). 
Dotted histogram is the original number-luminosity function of 
the volume limited sample of sources within d < 4.1 kpc. This 
function was multiplied by a correction factor 6.6, calculated from 
densities of HMXBs in different galactocentric annuli. Red solid 
histogram is a luminosity function of the whole sample (normal- 
ized to the number of HMXBs over the whole Galaxy), calculated 
taking into account luminosity dependent corrections (due to lim- 
ited sensitivity of the survey), dashed red curve — the best fit 
model of the luminosity function with parameters from Table 2. 
Hatched area shows the number-lu minosity function of all classes 
of HM XBs in our Galaxy from iGrimm. Gilfanov. fc Sunvaevl 
j2002h . 



our sample we have limited ourself by an axially symmetric 
distribution of HMXBs in the Galaxy. We have divided the 
Galaxy into annuli with radii R g < 2 kpc, 2 — 5 kpc, 5 — 8 
kpc, 8 — 11 kpc, 11 — 14 kpc from the center (the distance 
from the Sun to the Galactic center is expected to be 8.5 
kpc). 

We assume that: 

• the surface density of HMXBs (sre/kpe -2 ) is constant 
within each annulus; 

• shape of the luminosity function and its parameters are 
the same for all annuli. 

In our flux limited survey the sources with different lu- 
minosities can be detected within different galactic areas. 
In order to take this effect int o account we have adopted 
method 1/Vmax (|Schmidtlll968h . For estimation of parame- 
ters of the HMXBs lumi nosity func tion <j>{L) = dN/dL we 
used the Cash-statistics (|Cashlll979l ) as follows: 

C = 2 J2^J ^)5 , m ax J (L)dL-^ln[0(i l , J )5* max ,,(L l , J )] 

(1) 

Here, a summation j goes over the set of annuli and 



Parameter 



Value and la error 



10 36 erg s" 



1.40 ± 0.13(stat.)±0.06(syst.) 
> 2.2 

2.5± 2 J(stat.)±1.0(syst.) 



R g , kpc 



N{L > 10 35 erg s _1 ) kpc" 



0-2 




0.0±0.05(syst.) 


2-5 




0.1lt°; C O4(stat.)±0.02(syst.) 


5-8 




0.13l°;o3( stat -)±0.01(syst.) 


8-11 


(3.8 


t 2 {\) x 10~ 2 (stat.)±6.5 X 10" 3 (syst.) 


11-14 


(6.2 


t^a) x 10~ 3 (stat.)±4.8 X 10" 3 (syst.) 



a summation i goes over Nj sources within each annulus. 
Smax.j is the maximum area of annulus j, within which a 
source with the luminosity Lij can be detected. 

A minimization of the C-statistics gives us best fit pa- 
rameters of the luminosity function and its normalization 
in each annulus. We have adopted a simple broken power 
law shape of the luminosity function with slopes a\ and 02 
below and above the break at the luminosity L*: 



dN 
~dL 



Aj(L/L*y 



if L < L, 
ilL>L» 



(2) 



where Aj - normalization of the luminosity function in each 
annulus j. Best fit parameters of this model with uncertain- 
ties are presented in Table 2. Statistical uncertainties corre- 
spond to a la level; the systematical ones were calculated 
by variations of distances to sources as it was described in 
Section 3.1. 

The shape of the LF of wind-fed HMXBs demonstrates 
a break, or at least a curvature. From a purely statisti- 
cal point of view the statistical significance of the break 
is ~ 5%, or ~ 2cr (AC « 5.9 for 2 additional param- 
eters). Obviously, this can not be interpreted as a solid 
detection, but rather as a possible indication of its pres- 
ence. It should be noted additionally, that evidences for 
some fl attening of th e HMX Bs LF was mentioned ear- 
lier by IVoss fc Aiellol j201Ci ) for sources in our Galaxy 
and IShtvkovskiv fc Gilfanovl (|2005l ) for the population of 
HMXBs in the Small Magellanic Cloud. Thus it is likely 
that the gradual flattening of the HMXBs LF is real. 



3.3 Radial distribution 

The HMXBs surface densities (averaged over corresponding 
annuli) are presented in Table 2 and Fig. 5,6. It can be seen 
that the overall distribution of surface density of HMXBs in 
the Galaxy has a peak at galactocentric radii 2 — 8 kpc. 

Typical ages of high mass X-ray binaries are not 
longer then tens of Myrs, therefore it is natural to antic- 
ipate that their spatial density traces regions of the re- 
cent star form ation rather then the s tellar mass distribu- 
tion (see, e.g., lHaberl fc Sasa kl l2000l ). Previous measure- 
ments of the HMXBs distribution in our Galaxy (see e.g. 
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Figure 5. An illustrative view of the surface density of HMXBs 
in the Galaxy (the darker color of the annulus corresponds 
to the higher surface density of HMXBs). Black dotted and 
dashed curves show areas of the Galaxy, within which the IN- 
TEGRAL Galactic survey detects all sources with luminosities 
> lO 35,5 erg s — 1 and > 10 35 erg s -1 , respectively. Red dotted 
and dashed circles show distances, till which we detect all sources 
with luminosities higher than 10 35 erg s — 1 and 2 X 10 35 erg s — 1 
according to the adopted flux limit 0.7 mCrab in the inner part 
of the Galaxy —107° < I < 136°. Blue points indicate positions 
of HMXBs from our sample. All distances are in kpc. 



0.2 




Rodius, kpc 

Figure 6. Dependence of the HMXBs surface density (histogram, 
right axis) and star formation rate surface density (left axis) 
on the galactoccntric distance. Star formation rates are pre- 
sented by their upp e r and lo wer bounds (solid curves) from work s 
Guesten fc Mezgerl jl982h: iLvne Manchester, fc Tavloil <1985ft ; 



Chiappini. Matteucci. fe Romano! feOOlfT " Error bars on the his- 
togram represent statistical (larger error bars) and systematic 
(smaller error bars) uncertainties. Systematic uncertainties are 
estimated by variations of distances to sources. 



Grimm. Gilfanov. fc Sunvaevl 120021 : iLutovinov et alj 120051 . 
20071 ; iBodaghee et alj I2012T ) supports this general conclu- 
sion, but various incompleteness of previously available sam- 
ples of HMXBs precluded accurate estimates of their global 
density distribution. 

A comparison of the obtained HMXBs surface densities 
wi th the star formatio n surfa c e densities taken from papers 
of iGuesten fc Mezgerl dl982T); ILvne. Manchester, fc Taylor! 



\ 19851 ); IChiappini. Matteucci. fc Romano! (|200ll ) 
shows their very good correlation: N(HMXB, L x > 
10 35 erg s-^/kpc 2 « 5.5 x 10" 2 SFR/SFRq (see FigE]), 
here SFRq is the surface density of the star formation rate 
near the Sun. 

It is necessary to note, that on the average map, 
obtained with the INTEGRAL observatory, there are 5 
HMXBs with fluxes > 0.7 mCrab, but with unknown dis- 
tances (see Table [T]). These sources were not included into 
the calculation of the HMXBs surface densities and the lu- 
minosity function. This means that the true surface densities 
of HMXBs in the inner rings of the Galaxy (where most of 
sources are concentrated) might be slightly (~ 14%) higher 
than presented in Table 2. 

3.4 Vertical distribution 

We have fitted a vertical distribution of HMXBs in two com- 
monly used ways: 1) with a simple exponential model of their 
volume density p = poexp(— z/h) and 2) with a model of a 
self-gravitating isothermal disk p = po sech 2 (z / \^2h) . The 



former model gives the best fit height h = 85l 15 pc, the 
latter one h = 90 ± 15 pc. The systematic uncertainty of 
the scale-height due to a limited accuracy of distances to 
sources is smaller than the statistical one. This scale-height 
of the HMXBs distrib ution is somewhat smaller than tha t 
presented in pape r s of Grimm. Gilfanov. fc Sunvaevl (|2002h ; 
iDean et all (|2005l ) ; IBodaghee et all (|2007l ). likely due to the 
higher completeness and uniformity of our sample. 

It is important to note that the scale-height of the 
HMXBs distribution is larger than the one of the distri- 
bution of massive stars in the Gala xy, e.g., the sample o f 
WR stars has a scale of ~ 45 pc (IConti fc Vaccal 199 1 
OB star formation reg ions ~ 30 pc dBronfman et allpOQi 
open clusters ~ 50 pc (|Pandev. Bhatt. fc Mahrall98S ; Ijosh 



120051) . This indicates that HMXBs should have traveled a 
finite distance from their birth sites. If we will assume that 
HMXBs receive their systemic velocity during supernova ex- 
plosions, we can make a rough estimate of the kinematic age 
of HMXBs after the supernova explosion ( see, e.g., similar 
estimates in iBrandt fc Podsiadlowskllll995l ). 

Adopting the value ~ 50-90 kms -1 of the systemic ve- 
locity as a typical value for HMXBs (see e.g. iKaper et all 
1 19971 ; iHuthoff fc Kaperl |2002h we can estimate their kine- 
matic age, t — 50 pc/(50 - 90) kms -1 ~ 0.5 - 1 Myr. 
We should emphasize here that these values are applica- 
ble to the wind-fed population, which we study in our sam- 
ple. Roche lobe overflowing systems with luminosities above 
(0.5 - 1) x 10 37 erg/sec might have smaller ages. 
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4 PROPERTIES OF THE WIND-FED NS 
HMXB POPULATION 

Having collected the statistically clear sample of wind-fed 
binaries with accreting neutron stars, we now can try to 
understand physical parameters which determine their pop- 
ulation. 



4.1 Slice of HMXBs with fixed orbital periods 

The simplest picture of neutron stars acc reting from a stellar 
wind w as developed in classic a l works oflHovle fc Lvttletonl 
l|l94(t l; lBondi fc Hovlei (|l944r ); |Pavidson fc Ostrikerl (|l973h . 



In this framework the luminosity of the accreting neutron 
star simply depends on the mass flow, intercepted by its 
gravitational field. Masses of neutron stars lie in the narrow 
range of values (~ 1.3 — 1.9Mq), therefore they do not in- 
fluence strongly on the HMXBs accretion luminosity. The 
mass M, intercepted by a neutron star from the stellar wind 
of the companion mainly depends on three parameters: 1) 
the stellar wind mass loss rate M w , 2) the distance between 
components in the binary system a and 3) the stellar wind 
velocity v v (here we adopt that the wind velocity is much 
larger than the velocity of the orbital motion of binary stars; 
such an assumption is valid for the fast wind of young mas- 
sive stars of our sample). Roughly it can be expressed as 

M oc Af w ig 4 a~ 2 . (3) 

The orbital separation in the binary system depends on 
its orbital period. For a fixed distance between components 
in HMXBs the X-ray luminosity of the neutron star will 
be a function of the mass loss rate and the wind velocity 
of the optical star. Speaking very generally, for some fixed 
distance between companions in the binary system its X-ray 
luminosity (the accretion powered luminosity of the neutron 
star) should depend mainly on the mass of the normal star. 
The larger mass of the optical star should lead to the higher 
X-ray luminosity of the binary. Lets compare these simple 
arguments with properties of the real sample. 

In Fig[7]we present the cumulative luminosity function 
of the wind-fed HMXBs within a limited range of orbital 
periods from 4 to 10 days (the latter restricts somehow the 
distance between binary components). 

From this figure it is clearly seen three main properties 
of this sub-sample: 

• persistent HMXBs with orbital periods 4-10 days have 
luminosities in the range from 10 35 erg s _1 to 2 x 10 36 erg s _1 
and their luminosity function can be approximated by a 
power law dN/dL oc L~ 1,6 (dashed line); 

• there is a tentative trend that masses of optical stars 
(in cases where we can find them in the literature) decrease 
towards the lower X-ray luminosities; 

• the lowest luminosity of the HMXB in this small sample 
is L x ~ 0.8 x 10 35 erg s _1 . 

The power law shape of the X-ray luminosity fu nction of 
HMX Bs with fixed orbital periods was explained bv lPostnovl 
(|2003h by properties of the mass distribution of secondaries 
in the binary, and scalings of their stellar wind mass loss 
rate. 

An additional important hint from this plot is that for a 
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Figure 7. Cumulative luminosity function of wind-fed NS- 
HMXBs from our sample with orbital periods in the range of 
4 — 10 days. Dashed curve demonstrates a model prediction cal- 
culated from a simple analytical formula dN/dL oc L -1 - 6 in the 
luminosity range 2.2x 10 35 < L x < 8x 10 36 erg s _1 . The compan- 
ion star masses in some binaries are written along the histogram . 
They were taken fro m papers of Tomsick & Mutcrspaueh l|2010h 
for XTEJ18 5 5-026, IMason et al.l l l201lh for ICR J1802 7 -2016 , 
iLevine et all (120041) for X 1908+075 , iQuaintrell et alj 1120031 1 
for VelaX-1 jCox Kaper. fc Mokieml ll2005ri for 4U 1907+097. 
IMason et ail fcOlOn for EXO 1722-363. Dotted line shows an ap- 
proximate value of the X-ray luminosity of the HMXB with 
M 2 ~ 8 - 10M Q (see text for details). 



given orbital period, it seems that there should be a minimal 
X-ray luminosity of the wind-fed HMXB. 

This conclusion is not absolutely robust from the ob- 
servational point of view because at luminosities below 
L x ~ (2 — 3) x 10 35 erg s _1 our sample have a low complete- 
ness with the respect to orbital periods of binaries. Never- 
theless, usually there are no significant observational prob- 
lems to determine orbital periods in the range of 4-10 days, 
therefore we do not expect that this incompleteness is large. 

From the theoretical point of view this lower bound- 
ary of the X-ray luminosity of the wind-fed HMXB can be 
easily understood. Let's consider the accretion onto the non- 
magnetic neutron star, i.e. without an additional 'accretion 
flow-neutron star magnetosphere' interaction, which can in- 
hibit the accretion and thus significantly diminish the time 
average X-ray luminosity of such a binary system. If we will 
fix the orbital period of the binary system (that approxi- 
mately corresponds to the fixed distance between compan- 
ions), it will have the minimal X-ray luminosity if the opti- 
cal star has the lowest mass. We can expect that the lowest 
masses in HMXBs should be approximately 8 — IOMq. If 
we will use the 'X-ray lumin osity - mass of the optical com- 
panion' scaling from work of iPostn ov (2003) L x oc Af|' 5 , we 
can estimate the minimal X-ray luminosity of the HMXB in 
this small sub-sample. The resulted minimal X-ray luminos- 
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Figure 8. Orbital periods and mean X-ray luminosities of 
HMXBs from the flux limited sample (open circles). Gray filled 
circles denote positions of known supergiant fast X-ray transients, 
which still passed our selection criteria for the "persistency" due 
to their faintness. In addition to sources in the Milky Way we 
also demonstrate positions of two HMXBs in Magellanic Clouds, 
which accrete via the Roche lobe overflow of the optical star. Dot- 
ted horizontal line shows the level above which we know orbital 
periods for more than 90% of sources from our sample. Dashed 
line shows an approximate lower boundary of the "allowed" area 
for wind- fed accreting sources (see text). 



ity will be in a range of (0.7 — 1) x 10 35 erg s , that well 
agrees with the distribution, presented in Fig[7] 



4.2 Toy model of the wind-fed NS-HMXB 
population 

In this section we will try to simulate the global prop- 
erties of the population of NS-HM XBs, incorporating 
main ingredients from 



1940) 



1973) 



works 
(Il944tl; 



o f iHovle fe Lvt tlcton 



3ondi fc Hovld (ll944h; iDavidson fc Ostriker 
Lamers. van den Heuvel. fc Pettersonl 



Iben. Tutukov. fc Yungelson 
Bhadkamkar fc Ghoshl lj2012 ) 



(|l995T l; IPostnovl 



The main simplification in the modeling of the HMXBs 
population comes from the fact that their lifetimes are small. 
Thus, it can be securely adopted that properties of these 
systems depend only on their current parameters and are 
not significantly influenced by the previous history. In our 
simple toy model we will use this fact explicitly. 

As it was mentioned in Section 4.1 the X-ray luminosity 
of the neutron star fed by a fast stellar wind of a massive star 
is determined by: 1) the stellar wind density at the position 
of the compact object, and 2) the wind velocity. 

The mass loss rate in the stellar wind M w is mainly a 
function of the mass of the optical star M2 a nd the wind 
velocity w w (e.g. ICastor. Abbott, fc Kleinlll975h . 



(4) 



where L2 is the luminosity of the optical star, which by- 
turn is connected with its mass M2, e a dimensionless effi- 
ciency parameter. It depends on the type and temperature 
of th e optical star l|Lamers. van den Heuvel. fc Pettersonl 
Il976l) and can be varied in the ra nge of e ~ 0.4— 1.0 for typ- 
ical t emperatures of supergiants (|Chatv ; 2008; Rahoui ~et al.l 
120081 ). In following calculations we adopt e ~ 0.6. 

For a fast stellar wind the mass accretion rate on the 
neutron star can be estimated as follows: 



M « M v 



/ GM n 

y fw,ns* 



(5) 



where K w , ns - the wind veloci ty at the position of the neutron 
star v v , nB ~ v v ,(l-R 2 /a) 1/2 (|Castor. Abbott, fc Kl"em1ll975l : 
IVinkteOOOh , which might be significantly lower than the wind 
velocity at infinity u w if the neutron star is not far from the 
optical star, R2 is its radius. 

For the mass-radius relation for optical stars in 
the binary system we adopt R2/R0 ~ 0.9 AI2/ M©, 
which was calculated for sources in our sample. It 
approximately true in general (but with a large 
scatter due to evolutionary effects). For the mass- 
luminosity relation we adopt L2/LQ = 19(M2 /M ) 2 ' 76 
l|Vitrichenko. Nadvozhin. fc Razinkovall2007l '). 

Combining above formulas with the assumed mass to 
energy conversion during the neutron star accretion L x = 
O.lMc 2 , we can simply estimate that the X-ray luminosity 
of the wind-fed neutron star is: 



5.4fc x 10 d 



M 2 
10 Mi 



10 Rr. 



v v [1 — R2 /a] 
1000 km s" 



1/2 



(6) 



where coefficient k incorporates all uncertainties in the ap- 
proach used for the hard X-ray luminosity estimation. To 
describe the observed population of HMXBs (see Fig |8l9[) in 
the best way this coefficient should be around ~ 1.4 — 1.5. 
Taking into account a number of suggestions and simplifi- 
cations in our model it can be considered as a reasonable 
value. 

The line of reasoning, enlisted above, predicts that per- 
sistently accreting wind-fed neutron stars should occupy 
some particular region on a scatter plot P or b— L x (see Fig[8]). 
Indeed, at any given orbital period P or b (which can be trans- 
lated into the orbital separation) there should be a lower 
limit on the HMXB X-ray luminosity, corresponding (within 
the framework of our simple consideration) to the binary 
with the smallest masses of the companion star and thus to 
the lowest values of the mass loss rate. This lower bound 
should have a functional form L x oc a~ 2 oc P OI ^ 3 at large 
orbital periods and, roughly speaking, divide the P or b — L x 
diagram into two areas - "allowed" for wind- fed neutron 
stars and "forbidden" for them. It is remarkable that we do 
see this lower boundary on the diagram P or b — L x , shown on 
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Figure 9. Luminosity function of persistent HMXBs in the Milky 
Way as seen by INTEGRAL (crosses) along with the model of 
the simulated HMXBs population (solid line). Long-dashed and 
short-das hed lines show best fit models of the HMXBs LF ob- 
tained bv lVoss fc Aiellol fculCh and lGrimm. Gilfanov. fc Sunvaevl 
[120021 ), respectively. The former one was calculated in the 15 — 55 
keV energy band, where the luminosity of NS-HMXBs is practi- 
cally the same as the one on the 17 — 60 keV energy band. There- 
fore the corresponding LF is prese nted in its respective luminos- 
ity int erval. The LF in the paper of lGrimm. Gilfanov. fc Sunvaevl 
l|2002h was calculated in the 2 — 10 keV energy band. Therefore, to 
recalculate it to the 17 — 60 keV energy band we used an approxi- 
mate ratio between fluxes F2_io]j e v/Fi7-60keV — 0-5, which was 
derived from the analysis of b roadband spectra of t ypical HMXBs 
with neutron stars (see, e.g.. iFiliopova et al]|2005T l, 



FigE] (dashed line). This provides a significant support for 
our simple approach. 



4.3 Luminosity function of the simulated 
population of wind-fed HMXBs 

In order to simulate the population of HMXBs with neu- 
tron stars we need to assume some distribution o f masses 
of seco ndaries and orbital periods of these systems. iPostnovl 
suggested that secondary stars in HMXBs have an 
initial (Salpeter) mass distribution ignoring a possible evo- 
lution of this d istribution before the beginni ng of an active 
HMXB phase. iBhadkamkar fe Ghoshl (|2012T ) recently ana- 
lyzed this evolution and showed that really the overall power 
law slope of this distribution does not change significantly in 
the mass interval 15 — 4QMq. Therefore, for our subsequent 
modelling we adopt the mass distribution of secondaries in 
the form dN/dM oc M" 2 ' 35 in the mass interval 10-60 M@. 

For logarithm of orbital periods we assume a flat dis - 
tribution (see e.g. I6piklll924l ; iMasevich fc Tutukov|[l98Sl ). 
but restrict this distribution by the interval 1.3 < 
log P(days) < 3.3 with gaussian tails of the width 0.4 
in order to mimic results of the p re-HMXB evolution of 
binaries |Bhadkamkar fc Ghoshl lioia ). 



It is necessary to emphasize here that in our toy model 
we not take into account several effects, which might be 
important for the detailed comparison of properties of the 
simulated HMXBs population with the observed one. In par- 
ticular: 

• we adopt some specific mass-radius relation for optical 
stars; this assumption led to the fact that the wind velocity 
at the infinity w w is the same for all stars in the simulated 
population; an evolution of young massive stars can lead to 
a spread (at least) in this value; 

• we do not take into account that an X-ray illumina- 
tion of the wind matter from the accreting neutron star can 
lead to a deviation of the wind flow from a simple analytic 
formula (4); 

• we do not take into account that some part of the persis- 
tent NS-HMXBs population can accrete a matter from the 
dense equatorial wind of Be stars (i.e. XPer or X 1145-619); 
it will raise their X-ray luminosities above what might be 
expected for a binary system with a spherically symmetric 
stellar wind; 

• we do not take into account that depending on the mag- 
netic field and spin period of the accreting neutron star the 
accretion flow might be sto pped by a propeller mechanism 
l|lllarionov fc Sunvaev|[l975l ); it should lead to a disappear- 
ance of binaries from different parts of our P rb — £x diagram 
and can distort the shape of the luminosity function. 

Nevertheless, below we show that in spite of these limi- 
tations our toy model reasonably well describe observational 
appearances of HMXBs. 

A comparison of the luminosity function of the simu- 
lated wind-fed HMXBs population with that measured for 
such systems in our Galaxy is presented in FigO It is clearly 
seen that in spite of our simplistic approach to the simu- 
lation of the HMXBs population their luminosity function 
closely follows the observed LF of wind-fed NS-HMXBs in 
the Milky Way. At the same time, the measured luminos- 
ity fuiicHrjn_of_FIM2£Bs_Js^ 

of iGrimm, Gilfanov, fc Sunvaevl (|2002h and IVoss fc Aiellol 
(|2010t ) (dotted and dashed lines, respectively). The origin of 
this difference can be connected with the absence of black 
hole accretors, Roche-lobe filling systems and transients in 
our sample, and another way of the correction for the in- 
completeness. 

The main properties of the luminosity function of the 
HMXBs population are the steep cutoff at luminosities 
above ~ 10 erg s _1 and its flattening at luminosities 
below ~ iq34-34.5 gr g g -i According to our simple pop- 
ulation model it is likely that wind-fed HMXBs are lim- 
ited by a hard X-ray luminosity ~ 5 x 10 36 erg s _ . 
Thi s conclusion is not new and was ment i oned already by 
e.g. lLamers. van den Heuvel. fc Pettersonl (|l976l ). More lu- 
minous HMXBs should have a different nature. For exam- 
ple, they can harbor WR stars with more powerful winds, 
can accrete from Roche lobe overflowing massive giants 
(which can provide a mass accretion rate as high as 1O _4 M0 
yr _1 or more), can have black holes as primaries, etc. In 
our Galaxy we know only a few such sources (e.g., Cyg 
X-l, Cyg X-3), which, being taken together with wind- 
fed NS-HMXBs, will make the bright part of the HMXBs 
LF flatter than we see it on Fig(9] This fact could be a 
reason why the luminosity function of HMXBs in galax- 



12 Lutovinov et al. 



1000 



u 




QJ 




irt 




Ol 

L— 


100 


QJ 




in 








o 






10 


£ 












O 




c 


1 


E 




_5 




>v 




ra 


0.1 


1 

X 




T3 








Ha 


0.01 



SMCX-1 



MC X-4 . 

<S*V °A01657 GX301 . 2 



complete 
Porb 




0.001 



Orbital period, days 

Figure 10. Density of sources, produced with our toy model of 
the wind-fed NS-HMXBs population (gray scale) along with the 
distribution of known galactic sources on the plane P or b — L x 
(same as Fig[HJl. 



ies with high star formation rates continues with a power 
law d log N/d log L x ~ —1. 6 towards luminosities much 
higher than in our sample (|Grimm. Gilfanov. fc Sunvaevl 
120031 ; iMineo. Gilfanov. fc Sunvaevl I2012T ). The flat part of 
the LF at low luminosities (below ~ 10 34 erg s _1 ) is a con- 
sequence of a log-constant distribution of orbital periods of 
binary systems. But, it is necessary to note, that at the 
present level of the survey sensitivity this flattening is not 
very significant. 

The corresponding distribution of the simulated 
HMXBs population on the P or b — L x diagram is shown in 
Fig llOl by a gray scale. It is clearly seen that it forms the 
"allowed" region on the diagram, where most of persistent 
HMXBs are located. 



4.4 Supergiant fast X-ray transients 

Monitoring observations of the Galaxy with the INTEGRAL 
observatory have revealed a new pattern of a variability of 
more than dozen of high mass X-ray bin aries - short bright 
flares with long per i ods of quiescence JSmith et alj 19981; 
Sunvaev et all |2003t ISguera et all 120051 ; iNegueruela et al.l 



Sunyacy et al. ZUWJ; ogucra et au l/uuot UNcgucmcia et al.l 
20061 : ISidolill201lh . These sources rapidly have grown up to 



be a special subpopulation of HMXBs - supergiant fast X- 
ray transients (SFXTs, virtually all of them shown to be 
containing early type supergiant companions). 

An origin of such flares is not yet fully understood. 
Proposed models can be separated roughly in two main 
branches: a) flares occurring due to an occasional ac- 
cretion of blobs of the mat ter from the clumpy stellar 
wind of a supergiant (see, e.g., Walter fc Zurita Herasll2007l ; 



lOskinova. Feldmeier. fc Kretschmar|l2012l ). and b) magnetic 
arrest of the accretion due to a rotating neutron star 
magnetosphere and its occasional breakthrough (see, e.g., 
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Figure 11. Positions of supergiant fast X-ray transients (SFXTs) 
on the P or b — L x diagram. Here all sources are shown by two 
values - using the median and maximal fluxes, measured by the 
INTEGRAL observatory (solid circles). Dashed circle denotes the 
measurements of the mi nimal flux fo r IGR J16418-4532 based on 
the XMM-Newton data dSidolill20l3 ') and recalculated to the 17- 
60 kcV energy band. 



Grebenev & Sunvaev 


20071; Bozzo, Falanga. & Stellal 120081; 


Postnov et al. 20081; 


Shakura et al. 20121). Some combina- 



tion of these two types o f models was considered by 
iDucci. Sidoli. fc Paiziij (|2010l ). 

In order to reveal a true nature of the SFXTs phe- 
nomenon it is very informative to put these sources in a 
context of the outlined toy model of the wind-fed HMXBs 
population. 

A sub-sample of SFXTs (and SFXT candidates) with 
known distances and orbital periods is presented in Table [3] 
and shown in Fig llll Now every source has two values of its 
luminosity - the most probable (median) value, calculated 
using all INTEGRAL measurements, and the peak luminos- 
ity value. The last ones were extracted from corresponding 
light curves of sources binned into ~ 2 ksec bins. It should be 
noted here that the true peak luminosities might be higher 
in some cases if smaller time bins would be used. 

It is clear that the median luminosities of genuine 
SFXTs lie in the forbidden area, i.e. below the line of the 
minimal X-ray luminosity, allowed for the wind-fed NS in 
the binary system with the young star having the smallest 
mass loss rate in the wind (see Section 4.2). It explicitly 
means that most of the time the accretion in SFXTs is in- 
hibited by some mechanism. The only one source, which 
median value lies in the allowed area - IGR J16465-4507 - 
sometimes classified as an "intermediate" SFXT, and its flux 
variability is largel y related with the orbital modulation of 
the accretion flow (|Clark et al.ll2O10l ). Somewhat a similar 
situation is emerging for IGR J16418-4532, which is firmly 
established SFXT, but located near the border between al- 
lowed and forbidden areas. It happen only because we were 
limiting ourself with data of INTEGRAL/IBIS, which are 



Population of HMXBs in the Milky Way 13 



Table 3. Luminosities of some Supcrgiant Fast X-ray Transients.™ 



Name 


r b 


n X, peaky 


Distance, 


p 

r orb > 


References'* 




10 35 erg s _1 


10 35 erg s" 1 


kpc 


days 




IGRJ16418-4532 


9.8 


334.2 


13.0 


3.75 


1, 2 


IGRJ16465-4507 


1.4 


74.6 


9.4 


30.32 


3, 4 


IGRJ16479-4514 


0.31 


35.8 


2.8 


3.32 


3, 5 


IGRJ17391-3021 


0.04 


29.1 


2.7 


51.47 


3, 6 


IGRJ17544-2619 


0.09 


29.4 


3.6 


4.93 


7, 8 


SAXJ1818.6-1703 


0.024 


17.1 


2.1 


30.0 


9, 10 


IGRJ18483-0311 


0.22 


16.1 


2.8 


18.55 


9, 11 



in the 17—60 keV energy band 
most probable (median) luminosity 

maximum luminosity, averaged over ~ 2 ksec observations 

references correspond to the distance and orbital period measurements 



References: (D lllhaty elTa! f . ( 2) ll'or bet eTSQ W ). Ql lRahoul <,T5T\ mUSl . (4) |l'lark et ST} 1217111 ) ■ (5) IJain, Paul, fc Duttal t^OTSI ). (6) |t>raye ctST} |2T)Tg1 ). 
(7) Pdlim, Chaty. L Ne K „cruela ,2006). (8) Clark et al. ,2000). (0) Torrejon et al. ,2010), (10) Zurita Heras & Chaty ,2000). (11) Seuera e, al. ,2007). 



significantly worse that thos e of focusing X-ray telescopes. If 
we instead use the results o f lSidolil(|2012h . obtained with the 
XRT/Swift telescope, the median luminosity of the source 
in the quiescent state will be ~ 2 x 10 35 erg s _1 , that is ~ 5 
times lower than it was measured by INTEGRAL. In this 
case the source position in the quiescence state will be deeply 
in the forbidden area (dashed line and circle in Fig llll) . Fi- 
nally, concerning IGR J16418-4532, it is necessary to note a 
possible high uncertainties for its distance. 

The clumpy wind model can produce flares, but it can 
not explain why the median X-ray luminosity of the SFXT 
binary is much lower than that of the similar binary with 
the similar orbital period, but not classified as the SFXT. 

In order to see directly the outlined dichotomy we can 
compare sources IGR J17544-2619 and IGR J18027-2016. 
Optical stars in both binaries have similar spectral types 
- 091b for the former and OBsg for the latter, both have 
the similar orbital periods (around 5 days) and thus the sim- 
ilar separation between components in the system. But the 
first one is a fast transient (SFXT) and the latter one is a 
normal accreting pulsar. This example clearly shows that 
there are some additional parameters in the binary system 
(e.g., magnetic field and spin period of the neutron star), 
which determines its global behavior. 

The position of the SFXTs binary along the Y-axis in 
Fig llll should depend on the duty cycle of its activity. It is 
natural to propose that the duty cycle should be higher and 
the ratio of the peak fluxes to the median ones should be 
smaller if the source is closer to the allowed region. Such 
a source will "degenerate" into the normal persistent su- 
pergiant HMXB if its position will be within the allowed 
region. Note, that systematic studies of known SFXTs do 
show such a trend for all sources e xcept for IGR J16479-4514 
ijRomano et al.l l2009: Sidoli 2011). It is possible that such a 
deviation in the overall trend is related with the fact that in 
this binary the optical star is close to filling its Roche lobe 
and our toy model of the wind accretion is not applicable to 
it. 



5 FAINTEST SOURCES OF THE INTEGRAL 
SURVEY AND BEYOND 



A general understanding of the luminosity function of 
HMXBs and their Galaxy-wide distribution provides us a 
possibility to make predictions of a number of persistent 
sources at fainter fluxes. 

Using the real sensitivity of the INTEGRAL Galactic 
plane survey instead of our adopted value 0.7 mCrab (see 
Fig[2]) we can predict that in total we should detect up to 
~ 54 persistent HMXBs (here we have taken into account 
that surface densities, calculated in Section 3.2, do not in- 
clude 5 HMXBs with unknown distances and thus the true 
surface densities of HMXBs might be ~14% higher than it 
is presented in Table 2). Forty two of them were already 
observed with the adopted flux limit 0.7 mCrab. Below this 
limit there should be only ~12 persistent HMXBs among all 
unidentified INTEGRAL sources. Six of them were already 
detected in the survey (see Table [T]). Judging from estimates 
of surface densities of different types of sources on the sky we 
can anticipate that the majority of the remaining galactic 
sources will be cataclysmic variables (CVs). 

This conclusion follows from Fig |12l where the expected 
number of HMXBs in flux limited surveys of the Galacti 
plane is shown in line with ones for active galactic nuclei 
(AGNs) a nd CVs. A number-fl ux function of AGNs was 
taken from lKrivonos et alj (|2010t ). a number-flux function of 
cataclys mic variables was calculate d using the parametriza- 
tion from lRevnivtsev et al.l |2008bl ). It is clear that the num- 
ber of CVs and AGNs begin to dominate over HMXBs at 
fluxes below ~ 10~ 12 erg s" 1 cm" 2 , which will be achieved 
with t he new generation of hard X-ray tele scopes like NuS- 
TAR (|Harrison et al.ll2010T ) and Astro-H (|Takahashi et all 
2010). Small fields of view of these instruments will limit 
the detection rate of HMXBs. In order to increase the num- 
ber of known persistent HMXBs in the Galaxy (mainly due 
to low luminosity sources ) one need to use larg e survey mis- 
sions like Spectrum-RG (|Pavlinskv et al.ll2009h . 
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Figure 12. Surface density of HMXBs in the direction to the 
Norma (345° < I < 330°) and Scutum regions (15° < I < 35°) in 
the Galactic plane (gray areas). Lines show predictions of num- 
bers of different types of sources in these areas: AGNs (dashed 
line), CVs (dotted line) and HMXBs (solid lines). Predictions of 
a number of HMXBs at faint fluxes were done using two types of 
their LF with different slopes d log N/d log at luminosities below 
10 34 erg s — 1 (at higher luminosities the LF was taken as mea- 
sured, see Table 2): 1) an upper solid curve - the slope is fixed at 
the value —1.4, as measured at 10 34 erg s — x < L x < 10 36 erg s _1 
2) a lower solid curve - the slope of the LF is fixed at the value 
— 1, as predicted from our toy model (see Fig[9]l. 



6 SUMMARY 

We have constructed for the first time a well defined sample 
of non-transient High Mass X-ray Binaries in our Galaxy 
using the flux limited 9-years long survey of the INTEGRAL 
observatory. Our results can be summarized as follows: 

• The majority of persistent HMXBs accrete matter from 
the stellar wind of their supergiant /giant companions; 

• The luminosity function of wind-fed persistent HMXBs 
with accreting neutron stars can be described by a broken 
power law with the break around ~ 2 x 10 36 erg s _1 ; at high 
luminosities the power law slope of the differential luminos- 
ity function is 7bri g ht > 2.2, at low luminosities 7f a i nt w 1.4; 

• Using this luminosity function we have showed that the 
predicted number of HMXBs will be significantly lower than 
number of CVs and AGNs in future surveys with the better 
sensitivity (NuSTAR, Astro-H); 

• The spatial density distribution of wind-fed NS-HMXBs 
over the Galaxy have been measured. We have divided the 
Galaxy into several annuli and showed that the HMXBs sur- 
face density has a maximum on galactocentric distances 2-8 
kpc. Such a distribution correlates well with the distribu- 
tion of the surface density of the star formation rate in the 
Galaxy. 

• We have constructed a simple toy model of the wind- fed 
NS-HMXBs population and showed that it can adequately 
describe properties of the observed population. The model 
clearly shows that wind-fed HMXBs should disappear at lu- 



minosities higher than ~ 5x 10 36 erg s _1 and their LF should 
flatten at L x < 10 34 erg s _1 . The overall shape of the lumi- 
nosity function of the simulated population agrees with the 
measured one; 

• We argued that for wind-fed NS-HMXBs there is an 
"allowed" region on the P or b — £x diagram due to general 
properties of the wind accretion. We demonstrated that all 
persistent wind-fed HMXB lies in this allowed area. 

• All supergiant fast X-ray transients (SFXTs) in their 
quiescent state lie in the "forbidden" area of the P rb — L x 
diagram and only during flares they "jump" into the "al- 
lowed" region. This strongly supports the idea that their 
transient behavior is caused by some kind of (possible mag- 
netic) inhibition of the accretion onto the neutron star. 
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